Hyperbolic Chaos of Turing Patterns 
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We consider time evolution of Turing patterns in an extended system governed by an equation of 
the Swift-Hohenberg type, where due to an external periodic parameter modulation long-wave and 
short-wave patterns with length scales related as 1:3 emerge in succession. We show theoretically and 
demonstrate numerically that the spatial phases of the patterns, being observed stroboscopically, 
are governed by an expanding circle map, so that the corresponding chaos of Turing patterns is 
hyperbolic, associated with a strange attractor of the Smale- Williams solenoid type. This chaos is 
shown to be robust with respect to variations of parameters and boundary conditions. 
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In nonlinear dynamics the notion of structural stabil- 
ity, or robustness, is one of the key tools allowing one 
to specify systems and effects that are really significant 
for theoretical and numerical researches, and especially 
for practical applications P, 0| • Among chaotic attrac- 
tors, structural stability is intrinsic to those possessing 
the uniform hyperbolicity ("the systems with axiom A"), 
mathematical examples of which were advanced already 
since 60's - 70's [3j-|6j]. That time, such attractors were 
expected to be relevant for various physical situations 
(such as hydrodynamic turbulence), but later it became 
clear that the chaotic attractors, which normally occur in 
applications, do not relate to the class of structurally sta- 
ble ones. This is an obvious contradiction to the principle 
of significance of the robust systems mentioned above. 

Recently, this inconsistency has been partially resolved 
by introducing a number of physically realizable systems 
with hyperbolic chaotic attractors [7HlO|. It has been 
shown that simple systems of coupled oscillators that 
are excited alternately (in time) possess hyperbolic at- 
tractors of Smale- Williams type (for experimental real- 
izations, see )• Hyperbolic chaos in these systems 
is related to the dynamics of the phases of the oscillators, 
evolution of which on the successive stages of activity is 
governed by a Bernoulli-type expanding circle map. 

In this letter we develop a similar approach, but deal 
with the spatial phases of patterns in a spatially extended 
system. We demonstrate the occurrence of hyperbolic 
chaos in dynamics resulting from an interplay of two Tur- 
ing patterns of different wave lengths arising in succes- 
sion. This advance, first, extends a toolbox for design 
of models manifesting robust chaos. Second, it suggests 
a novel direction for search of situations associated with 
hyperbolic chaos in the context e.g. of fluid turbulence, 
convection, and reaction-diffusion systems. Third, the 
description in terms of truncated equations for ampli- 



tudes of spatial modes leads to new prototypical low- 
dimensional model systems with hyperbolic attractors. 
(Note analogy with the Lorenz equations, which were 
derived originally as a finite-dimensional model for fluid 
convection.) 

Let us illustrate the approach with a concrete exam- 
ple based on the one-dimensional Swift-Hohenberg equa- 
tion (l2^ . Consider its following modification: 



d t u + [1 + n 2 (t)d 2 x ] 2 u = \A + Bx(x)]u - u 3 



(1) 



Here A is a positive parameter controlling the Turing 
instability. An instant value of k determines the wave 
number of the unstable Turing mode. In our case n(t) is 
assumed to be a periodic function: n(t) = 1 for nT < t < 
(n+l/2)T, and n(t) = 1/3 for (n + l/2)T < t < (n+l)T. 
This switching provides the excitation of two distinct al- 
ternating in time Turing patterns characterized by the 
dominating wave numbers, k = 1 and k — 3, respec- 
tively. The time interval T between the switchings is 
supposed to exceed the characteristic time duration of 
the formation or decay of the Turing patterns. A nonlin- 
ear cubic term in the equation is responsible for satura- 
tion of the instability. Also, the coefficient at the linear 
term in the equation is assumed to depend on the spatial 
coordinate that corresponds to the presence of a spatial 
non- uniformity characterized by a function x( x )'i its role 
will be clarified below. Assuming the ring geometry and 
periodic boundary conditions u(x, t) = u(x+L, t) (PBC), 
it is natural to set the length of the system as L = 2irl, 
with integer £, to get the geometry supporting the Turing 
patterns of both the wave numbers k — 1 and 3. 

The system operates as follows. In each time interval 
with Ac(t) = 1 the Turing pattern with the dominating 
wave number k = 1 arises, which is characterized by some 
spatial phase ip: u ~ U\ cos(x + (p) + U3 cos(3x + 3(f), 
where U3 <C U±, and U± is of the order of \f~A~. (The 
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Figure 1. (Color online) Hyperbolically chaotic patterns in 
model © at A = 0.6, B = 0.03, T = 25, L = 4tt, N = 64, 
PBC. 



third harmonic appears naturally due to the cubic non- 
linear term in the equation.) After the switch to k = 1/3 
the system becomes unstable in respect to the harmonic 
component with k = 3, while that with k = 1 starts 
to decay. The initial stimulation of the short-wave pat- 
tern is provided by the component C/3; so, it accepts the 
spatial phase 3<p. At the end of the considered time in- 
terval the first harmonic component practically disap- 
pears, and we have u ~ U^cosi^x + Stp), with U3 of 
the order of \/~A. After the next switch, when n = 1 
again, the third harmonic decays, but the first harmonic 
becomes unstable and starts to grow. A germ for this 
growth is provided by a component at the wave num- 
ber k — 1 arising from the combination of the decaying 
short-wave pattern and the spatially dependent coeffi- 
cient x( x )- If the Fourier expansion of x{x) contains a 
dominating second harmonic k = 2, the long-wave mode 
will arise with the phase 3<p, due to the term proportional 
to cos 2x cos(3.t + 3(p) = (1/2) cos(.t + 3<p) + . . .. Thus, on 
each complete period of modulation T the phase of the 
spatial pattern undergoes the tripling (up to a constant 
phase shift): <p n +i = 3(p n + const. This is an expand- 
ing circle map with chaotic behavior characterize d by 
the positive Lyapunov exponent A = In 3 ps 1.0986 \13\ . 
Since the phase map is uniformly expanding, the strobo- 
scope map corresponding to the transformation of the 
states u n (x) = u(x,t n ) from t n = nT + const to t n +\ is 
expected to be hyperbolic. 

Of course, this mode of operation occurs under the 
proper choice of the parameters. A value of A is selected 
to get an instability at k — 1 with a decay at k = 3, or 
vice versa, at successive half-periods of parameter modu- 
lation. The term B\(x) must be small (comparing to the 
fully developed pattern amplitude) to contribute only as 
a germ for the formation of the long-wave pattern, al- 
though this germ should be of a sufficient level to start 
the process with saturation on the time scale T. In fact, 
the requirements are not very strict: the described type 
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Figure 2. (Color online) Diagrams for spatial phases of Turing 
patterns at t n = (n + 1/4)T for A = 0.6, B = 0.03, T = 25. 
(a) Numerical solution of the system fTJ, PBC, L = 4tt, N — 
64. (b) Amplitude equations ©. (c) System {U with ZBC, 
L = 87r, N = 128. 

of behavior occurs in a fairly wide parameter range. 

Figure[T]illustrates the spatio-temporal behavior of the 
system observed for the case of PBC. The 3D-plot u(x, t) 
is obtained using computations on a spatial grid with 
the node separation Ax- = L/N, where N is a number 
of the nodes. One can observe the alternating evolu- 
tion of the Turing structures: a long-wave pattern first 
appears, then decays, and is replaced by a short-wave 
one. After the period T, the long-wave pattern reappears 
but with a different spatial phase (shift along x-axis), 
and the process repeats. As we show in Fig. Uta), the 
spatial phases recorded stroboscopically follow a chaotic 
map of the expected type. To obtain this diagram, we 
determine the spatial phases at t„ = (n + 1/4) T as 
tp n = &Tg[u(L/2,t n ) + id x u(L/2,t n )], where the spatial 
derivative d x u is estimated by the numerical differenti- 
ation, and the results are plotted in coordinates tp n +i 
versus (p n . This empirical map is of the expected topolog- 
ical type: one revolution for the pre-image corresponds 
to three revolutions for the image. 

To characterize chaos quantitatively and demonstrate 
its robustness, we calculate the Lyapunov exponents. 
Figure [3ja) shows the first five Lyapunov exponents for 
the stroboscopic map as functions of the parameter A. 
The chaotic mode of operation occurs above some thresh- 
old around A as 0.38. In the chaotic regime there 
is one positive Lyapunov exponent, which remains al- 
most constant in a wide parameter range. In partic- 
ular, at A = 0.6 the Lyapunov exponents are A = 
{1.018, -9.34, -9.34, -11.42,-18.64, . . .}. As expected, 
the largest exponent is close to In 3. As seen from the 
diagram, all the exponents depend on the parameter 
smoothly, without sharp spikes or dips. This is a mani- 
festation of robustness of the hyperbolic chaos 0, & E3| ■ 
The Kaplan- Yorke dimension of the attractor varies 
slightly, see the solid line in Fig. [3Jb); in particular 
D KY f« 1.11 at A = 0.6. 

To confirm the validity of the used spatial discretiza- 
tion, in Figure [3]Jc) we show the sixteen largest Lya- 
punov exponents obtained at a fixed length L with dif- 
ferent sizes of the numerical mesh N. The decrease of 
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Figure 3. (Color online) (a) Five largest Lyapunov exponents 
vs. A for the stroboscopic map of the system ([T]) at t n = 
(n + l/4)T, PBC. (b) Kaplan- Yorke dimension for PBC (solid 
line) and ZBC (dotted line), (c) PBC, first sixteen exponents 
for different N: pluses, crosses, stars and squares refer to 
N = 64, 128, 256, and 512, respectively, (d) ZBC, first five 
exponents. Other parameters are the same as in Fig. Q] for 
PBC and B = 0.03, T = 25, L = 8tt, N = 128 for ZBC. 



Ax = L/N corresponds obviously to approaching the 
continuous limit. The left-hand parts of the curves over- 
lap perfectly; so, the larger exponents are in good cor- 
respondence for all tested step sizes. The discrepancy 
visible in the right-hand part of the plot for large nega- 
tive exponents decreases with the growth of N. Hence, 
we can be sure that the properties revealed in the com- 
putations with the finite discretization size are valid for 
the continues system as well. 

Next, we perform a direct test of the hyperbolicity. 
The hyperbolicity implies that there are no tangencies be- 
tween the stable and unstable manifolds of orbits belong- 
ing to the attractor. Occurrence of a tangency is deter- 
mined by the zero angle between the expanding and con- 
tracting tangent subspaces spanned by the corresponding 
covariant Lyapunov vectors [l4j |. Following the method 
for testing hyperbolicity described in (l5| . we examine 
the distribution of these angles by considering the or- 
thogonal complement to the contracting subspace, which 



is normally much less dimensional then the contracting 
subspace itself. If there are K expanding directions, it 
is sufficient to calculate K orthogonal backward and for- 
ward Lyapunov vectors, to construct a K x K matrix P 
of their scalar products, and to check how close to zero 
is the normalized characteristic number 



d K = |det(P)|. 



(2) 



By the definition, < dx < 1. The procedure is applied 
at a representative set of points on a trajectory on the at- 
tractor. The distribution of cLk separated well from zero 
means that the chaos is detected as hyperbolic: the tested 
trajectory does not contain any points with tangencies of 
the expanding and conr acting Lyapunov vectors. 

In application to the stroboscopic map of the sys- 
tem (JTJ the calculations are simple because K = 1. For 
the parameters used in Fig. Q] we processed 10 5 points 
and observed that (1 — 5 x 10~ 5 ) < d\ < 1. It means 
that the expanding direction is always almost orthogonal 
to the contracting subspace. Thus, the conjecture that 
the attractor is uniformly hyperbolic is confirmed, but, 
of course, a rigorous mathematical proof of the hyperbol- 
icity would be desirable anyway. 

As in the system only two modes with the wave num- 
bers k = 1 and k = 3 are basically involved, one 
can expect that the essential properties of the dynam- 
ics can be described with a truncated model. To de- 
rive the low-dimensional model we proceed as follows. 
Accounting for the relevant modes, we use the ansatz 
u = a±(t) cosx+bi(t) svax+a,3(t) cos 3a;+& 3 (i) sin 3a;. As- 
suming x( x ) — cos 2a;, after the substituting to Eq. (JTJ, 
we multiply the resulting expression by cos a; and sin a:, 
and by cos 3a: and sin3x, and for each case perform the 
integration over the spatial period 27r. The result is a set 
of equations for the amplitudes of the modes, which can 
be compactly expressed in the complex form as 

ci =MiCi - i [3(| Cl | 2 + 2|c 3 | 2 ) Cl - 2Bc 3 + {Zc\c 3 - 2B)c\] 



c 3 =M3C 3 - 7[3(|c 3 | 2 + 2| Cl | 2 )c 3 - 2Bc x + c 



lJ • 



(3) 



where the asterisk denotes complex conjugation, c\ = 
oi + i&i, c 3 = a 3 + i& 3 , Hi = A - (1 - k 2 ) 2 , /z 3 
1 : 1 9k 2 ) 2 , and k = n{t). as before. Notice that the 
structure of the equations resembles that for the ampli- 
tude equations obtained for other models with hyperbolic 
attractors of Smale- Williams type d, S O • 

Figure |4] illustrates the dynamics of the model ((3]). 
Observe the switchings after each next half-period T/2. 
Different heights of the humps of Re c\ and Im c\ arise 
due to the variations of the phases of c\. The phases 
transform stroboscopically according to the triple ex- 
panding circle map; see Fig. Efb) for the diagram for 
the phases computed as = arg[ci (£„)]. The Lya- 
punov exponents evaluated for the stroboscopic map of 
the model © at A = 0.6, B = 0.03, T — 25 are 
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Figure 4. (Color online) Solution of Eq. J3]) at A = 0.6, B = 
0.03, T = 25. Solid (red) and dotted (blue) lines refer to 
ai = Re ci , and b\ = Im ci , respectively. 



A = {1.083,-12.5,-804.7,-806.5}, and the Kaplan- 
Yorke dimension is 1.09. Notice that the first expo- 
nent is close to In 3. The hyperbolicity test described 
above again shows that the expanding and contracting 
subspaces are almost perfectly orthogonal. 

Now we intend to demonstrate that the hyperbolic 
chaos can be observed in geometries distinct from the ring 
one. In particular, let us examine simple zero boundary 
conditions (ZBC): u(x, t) = for x < and x > L, which 
could, in principle, block the mechanism of the chaotic 
transformation of the phases (because near the ends the 
spatial phase is dominated by the boundary conditions). 
However, computations show that such a blocking occurs 
only in short systems. If the length is large enough, pat- 
terns in the middle part of the system still interact in the 
same way as for PBC, while the parts close to the ends 
undergo deformations to fit ZBC; see Fig. [2jc) and [3jd) 
for L = 8tt. The map for the spatial phases at the middle 
part of the system agrees well with the expected form. 
Moreover, there is a large parameter interval, where the 
system has a single positive Lyapunov exponent of value 
almost independent on A. At A — 0.6 the Lyapunov ex- 
ponents are A = {1.047, -1.59, -3.92, -4.97, -6.16, . . .}, 
and the Kaplan- Yorke dimension is 1.66. The hyperbol- 
icity test shows pronounced separation of d\ from the 
origin, although the distribution is wider (0.93 < di < 1) 
than for PBC. 

Summarizing, in this letter we have shown how the 
hyperbolic chaotic dynamics can emerge in extended sys- 
tems due to an interplay of spatial patterns with different 
wave lengths. In our model system the spatial phases of 
the patterns evolve in time according to the Bernoulli- 
type tripling map, and their dynamics is strongly and 
robustly chaotic, while the amplitudes behave in a rather 
regular manner. The mechanism of the hyperbolic chaos 
is similar to that in alternately excited oscillations, stud- 
ied earlier Q. In some respects, the chaotization of spa- 
tial phases appears to be easier for implementation (there 
is no necessity to have more than one involved subsys- 
tem). We have demonstrated the expected chaotic be- 
havior in the partial differential equation of the Swift- 
Hohenberg type, and in the truncated model represented 
by a set of ordinary differential equations. It should 



be emphasized that the kind of dynamics we consider 
is not specific for the Swift-Hohenberg equation only. 
Ingredients needed for the phase multiplication mecha- 
nism, namely, the alternation of patterns due to param- 
eter modulation, the non-linearity, and the spatial inho- 
mogeneity can be either found or created in many spa- 
tially extended systems. As expected, these results open 
prospects for the search and constructing for hyperbolic 
chaos in pattern-formation for systems in fluid dynam- 
ics (Faraday ripples, convection rolls) and in reaction- 
diffusion systems (Turing structures, advection induced 
patterns) [l2j]. In the case of microfluidic systems (l7l[l8l ] . 
an interesting question for future studies is the effect of 
hyperbolic chaos on Lagrangian mixing properties. 
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